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ABSTRACT 

In this paper, the stability of the time dependent differential equation for HIV infection of CD4+ T cell model 
was studied using the fractional order derivatives. It is discovered that the eigenvalues showed that the 
system was unstable with three different equilibrium points using Matlab 7.5.0.342 R2007b version. 


1. INTRODUCTION 
Weli et.al. (2014) studied the time dependent ordinary differential equation model on HIV 
infection of CD4+ T cells of a patient. The following three differential equations describe the 


model. 
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T(t)is the number of target cells, I(t) is the number of infected cells, V(r) is the viral load of 


the virions and the R is positive constant. œ is the maximum proliferation rate of the target 


cell, Gis the infection rate constant, dis the death rate of cells infected, pis the 


reproductive rate of infected cells, sis the rate at which the new T cells are created from 


source within the body, c is clearance rate constant of virions, d is the death rate of T(t) 


cells of CD4+ T cells and T x is the T population density at which proliferation shuts off. 


( -7 describes the logistic growth of the healthy CD4+T cells. 


max 
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2. LITERATURE REVIEW 
In this section, we shall review some existing literature done by some researchers. 
Arafa et.al. (2016) studied the stability of equilibrium points of a model of HIV infection of 
CD4+ T using Mittag-Leffler method. The results showed that the fractional ordinary 
differential equations has more advantages than classical integer-order modeling. 
Arafa et.al. (2013) studied the nonlinear differential equation system such as Lorenz 
system of fractional order. They used generalized Mittag-Leffler to solve the approximate 
and analytical solutions of the system. The result showed that the fractional order 
differential equation spread at a faster rate than the classical differential equation. 
Arafa et.al. (2012) studied the effect of changing the average number of viral particles with 
distinct sets of initial condition. They used the Generalized Euler Method. The result 
showed the numerical solution of the HIV primary infection. 
Goreno and Mainardi (2000) reviewed a variety of fractional evolution processes. They 
used fractional order whose solutions turned out to be related to Mittag-Leffler-type. 
Hegazi et.al (2011) studied the synchronization of chaotic and hyperchaotic systems using 
fractional order. The result showed the effect of the method used on the synchronized 
systems. 
Matignon (1996) studied the internal and external stabilities for finite-dimensional linear 
fractional differential systems. He used a fractional differential system in polynomial 
representation to establish the stabilities of the roots of some polynomial which lie outside 
the closed angular sector. 
Merdan (2007) studied HIV infection model of CD4* T cells using homotype perturbation 
method. His result showed the reliability and simplicity of the method. 
Merdan et.a/ (2011) studied HIV infection model of CD4* T cells using Variational Iteration 
Method. Their results showed the reliability and simplicity of the method. 
Rida and Arafa (2011) studied a linear fractional differential equations using the Mittag- 


Leffler function method. They results showed that method was efficient. 


3. FRACTIONAL ORDER DERIVATIVES 
Now, following the Arafa et.al. (2016), we introduce the fractional order derivative of HIV 
infection of CD 4+ T cells: 


D” a=s-ar san1- —)- rv 


max 


D® (1) = PTV -A (3.1) 
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D” (V)= pl -cV 
Where 77,,77,,77; > Owith initial conditions T (0) =0.1, 7(0) =0, V(O) =0.1. 


4. EQUILIBRIUM POINTS 
We suppose that: 
TtHh=x,) 
I(t) = x,(t) (4.1) 
V(t) = x(t) 
Thus for the equilibrium point, set D’x,(t) =0 for i =1,2,3.... 


Hence, (1.1) becomes: 


s—dx, van[ ZL }- Bx, x, =0 (4.2) 
Px,x, — Ox, =0 (4.3) 
PX, — cx, =0 (4.4) 


Solving (4.2), (4.3) and (4.4) to have that if x, =0 > x, =0. 





af, -dT x +d°T2, —2daT2, + a°T2, +4asT.,,, 
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Then, 
E OL T ol Ae doh 
(xe, x2, xi?) = ,0,0 
2a 
(4.5a) 
and, 
big g\ | Chg sal P12 <2dar 40 T?, + AasT.. 
(xe x% x)= ,0,0 
2a 
(4.5b) 


are equilibrium points. 


Next suppose x, # 0 and x, #0, then another equilibrium point is 





ed 


(cet, x x") = c spp, — dco, + ACT. — CS” SFT yy, — ACT yy, + LCT yy — = 
pp BOOT rax BCT as 
(4.6) 
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5. NUMERICAL STABILITY OF THE SYSTEM 


With the given parameter values: 
a = 6.8day™, B =0.0002mm~, 5 =0.5day', p = 1000mm day", 
s =(5day')(mm™~),c =10day",d =0.0lday",T.,,. = 1300mm day”. 


So, from (4.5a): x =1.2988e +003. 
Next, to form the Jacobian matrix, let F(x,,x,,x,), G(x,,x,,x,;)and H(x,,x,,x,)be three 


continuous and differentiable functions such that: 


F(x,.x).x,) = 9—dx, vas 2) (5.1) 
G(x, X3, X3) = AXX — X, (5.2) 
H(X, X3, X3) = PX, — CX; (5.3) 
So that the Jacobian matrix becomes: 
— 6.7977 0 —0.2598 
A= 0 —0.500 0.2598 
0 1000 —10 


Its eigenvalues are A, =-6.7977 <0, A, =11.5536>0, A, = -22.0536 <0. 


Hence, the equilibrium point 








oo ea cea ( ug, — Wg, +f T —2da??,, + T2, + 40ST... 
(xe, x27, x 1) = 


3 oo) is unstable. 
2a 








aT ax ~ AT, yd T, —2daT~., tA T + 4asT. 
Next, for: eraras) max max J max a max max max ,0,0 
a 


Then the Jacobian matrix becomes: 


6.7977 0 1.4719e —004 
A=| 0 —0.5000 —1.4719e—004 
0 1000 —10 


Its eigenvalues are 1, = 6.7977 > 0, A, =-0.5155 <0, A, =—9.9845 <0. 


Hence, the equilibrium point 








E OD e — AT yy, — fT, 2d, KAT + 40ST ye, 


,0,0 | is unstable. 
2a 
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Next for the point: 


erara SPBT vox ~ ACT roy + LCT yo -acs SPBT vex ~ ACT y + CT mon a 
pp BOPT rax BOOT rax 





The Jacobian matrix becomes: 


6.4588 0 —0.0050 
A=| 0.0696 -0.5000 0.0050 
0 1000 —10 


Its eigenvalues are 1, = 6.4556 > 0, A, =0.0051>0, A, = —10.5020 <0. 

Hence, the equilibrium point 

( nD aed e) z SPPT.,. — dcoT.,,. + ACT ax —a°C°S’ SPBT ax — dcoT yx + ACT ax =e) 
x; 1X5 > X3 = > > 


max 


pp BOLT rax BOLT rax 





is unstable. 


6. CONCLUSION 


It was studied in this work that the given HIV infection of CD4+ T cells model was unstable 


with the help of eigenvalues given the values of the parameters a, P, ô, p, s,c,d, T x- 


For a unique equilibrium point, the system gave at least a positive eigenvalue. There is a 
need for varying the data parameters or looking for the second order fractional derivatives 


of the discussed model. 
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